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The small-scale dynamo is a process by which turbulent kinetic energy is converted into magnetic 
energy, and thus it is expected to depend crucially on the nature of the turbulence. In this paper, 
we present a model for the small-scale dynamo that takes into account the slope of the turbulent 
velocity spectrum v(£) oc £®, where £ and v(£) are the size of a turbulent fluctuation and the typical 
velocity on that scale. The time evolution of the fluctuation component of the magnetic field, i.e., the 
small-scale field, is described by the Kazantsev equation. We solve this linear differential equation 
for its eigenvalues with the quantum-mechanical WKB-approximation. The validity of this method 
is estimated as a function of the magnetic Prandtl number Pm. We calculate the minimal magnetic 
Reynolds number for dynamo action, Rm cr it, using our model of the turbulent velocity correlation 
function. For Kolmogorov turbulence (i? = 1/3), we find that the critical magnetic Reynolds number 
is Rm^ it w 110 and for Burgers turbulence (i? = 1/2) Rm^ rit w 2700. Furthermore, we derive that 
the growth rate of the small-scale magnetic field for a general type of turbulence is P oc Re' l_1 '^ l+ '^ 
in the limit of infinite magnetic Prandtl numbers. For decreasing magnetic Prandtl number (down 
to Pm > 10), the growth rate of the small-scale dynamo decreases. The details of this drop depend 
on the WKB-approximation, which becomes invalid for a magnetic Prandtl number of about unity. 

PACS numbers: 47.65.Md,95.30.Qd,47.27.-i,84.60.Lw 



I. INTRODUCTION 

Magnetic fields are observed in the whole Universe on 
different length scales and they play an important role 
for virtually all astrophysical objects. To probe the ori- 
gin of these fields, it is crucial to gain better theoretical 
insight into their dynamical properties. 
Observations suggest that magnetic fields were already 
space filling in the early Universe. Hints toward high- 
rcdshift magnetic fields come from the study of Bernet 
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et al. pp. Recent studies using 7-ray observations of the 
intergalactic medium find a lower limit of Bq ss 10 — 
1CU 6 nG gHl], indicative of an early generation scenario. 
For the very early Universe, we have only upper limits, 
but no direct observations showing the presence of mag- 
netic fields. Yamazaki et al. [5] derive an upper limit of 
the magnetic field strength from the cosmic microwave 
background temperature anisotropy. They predict that 
the magnetic field strength at the present scale of 1 Mpc 
[B] is Bq < 4.7 nG. Moreover, Schleicher and Miniati [7 
present a method to derive the magnetic field strength 
in the reionization epoch (redshift z > 6). They show 
that the upper halo mass is set by the magnetic pressure 
and conclude, by using data from the reionization epoch, 
that the co-moving magnetic field strength is Bo < 3 nG 
(see also [8]). 

Besides these observational constraints, there are differ- 
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ent theories that describe the origin of these primordial 
magnetic fields. The first seed fields could already have 
been produced during inflation. Turner and Widrow [D] 
find that B Q ps fO _25 -fO _1 nG on a scale of 1 Mpc can 
be produced when the conformal invariance is broken. 
Following Sigl et al. [10], there is also the possibility of 
creating a magnetic field during the phase transitions in 
the very early Universe. They predict a field strength 
Bo w 10~ 20 nG from the electroweak phase transition 
and Bo w 10 -11 nG from the quantum chromodynamics 
(QCD) phase transition. 

The observed magnetic field strengths in high-redshift as 
well as in present-day galaxies are, however, typically or- 
ders of magnitude higher than predicted by the theories. 
For example the Milky Way has locally a magnetic field 
of (6 ± 2) x 10 3 nG [11] . Hence there must be additional 
and strong amplification. The simplest way of amplify- 
ing a magnetic field by gravitational compression under 
the condition of flux freezing, when the magnetic field is 
perfectly coupled to the gas. With the conservation of 
magnetic flux, one can directly show that B cx p 2 / 3 for 
spherical collapse, where p is the mass density of the gas. 
However, magnetohydrodynamical dynamos can amplify 
small initial magnetic fields exponentially in time, lead- 
ing to a growth of B much faster than /? 2 / 3 [12TH4] . A 
dynamo, in general, is a process by which kinetic energy 
is converted into magnetic energy. One has to distinguish 
two different types of dynamo: the large-scale dynamo, 
which is excited by large-scale motions, and the small- 
scale dynamo, which is excited by turbulence on very 
small scales. Dynamos amplify the magnetic energy ex- 
ponentially. The small-scale dynamo is important in par- 
ticular during the formation of the first stars and galax- 
ies, where turbulence is released from the gravitational 
potential well, providing the means to amplify magnetic 
fields [14]. The same processes may also occur during 
the formation of large-scale structure [15]. In most as- 
trophysical objects, the growth rate of the small-scale 
dynamo is very high and that can explain the high mag- 
netic field strengths observed in the primordial as well as 
in the present-day Universe. 

The small-scale dynamo was first suggested by Batchclor 
in 1950 [IB]. In 1967, Kazantsev developed a spectral the- 
ory of this process [17j[IB]- He derived a differential equa- 
tion, the Kazantsev equation, which describes the evolu- 
tion of the magnetic energy of the small-scale field. This 
equation has been solved for different special cases. For 
example, Subramanian |19j used this theory to describe 
the small-scale dynamo in the case of incompressible tur- 
bulence, i.e., Kolmogorov turbulence. In the limit of infi- 
nite magnetic Prandtl number Pm, which is the ratio of 
the magnetic to the hydrodynamical Reynolds numbers, 
Rm and Re, he found a growth rate of the small-scale 
magnetic energy of (15U)/(4L)-\/Re [2D]. Schekochihin 
et al. [21] analyzed the Kazantsev theory for an arbi- 
trary compressibility (and spatial dimension), also in the 
limit of large Prandtl number. However, they did not 
give a model for the correlation function of the turbu- 



lent velocity field, but used a small-scale approximation 
of this tensor. There are also alternative approaches for 
the limit of small Prandtl number based on the separa- 
tion of scales [22H24"] , which find similar solutions of the 
Kazantsev equation. Further ways to solve this equation 
exist, but go beyond the scope of this work (see for ex- 
ample [25] HBJ). 

In this paper, we aim to explore dynamo amplification for 
turbulence models with v(£) cx £ in the inertial range, 
i.e., the range of the kinetic energy cascade. Here, I is 
the size of a turbulent fluctuation and v{£) the typical ve- 
locity on this scale. The range of d goes from 1/3, which 
describes the incompressible Kolmogorov turbulence [27] , 
to 1/2 for Burgers turbulence [28]. In astrophysical ob- 
jects we expect that the exponent $ lies between those 
extrema 29-33 . Burgers turbulence is often referred to 
as highly compressible turbulence, obtained in the limit 
of very high Mach numbers. 

Federrath et al. [33] explore the Mach number depen- 
dence of the small-scale dynamo and show, that the dy- 
namo works for both, subsonic and highly supersonic 
turbulence. Furthermore, these numerical experiments 
indicate that the small-scale magnetic field saturates af- 
ter a certain time. They find that the saturation levels 
and growth rates depend strongly on the nature of the 
turbulence. However, most simulations typically have a 
magnetic Prandtl number of Pm ps 1, requiring further 
efforts and understanding for Pm 3> 1 and Pm <C 1. 
Here, we present an improved analytical model for the 
small-scale dynamo that takes into account highly com- 
pressible turbulence also. In Sec. [IT] we outline the main 
ideas and equations behind the Kazantsev theory and de- 
scribe the small-scale dynamo mathematically. We use 
the WKB approximation as a solution of the Kazant- 
sev equation, which is valid in the limit of large mag- 
netic Prandtl number. However, the evolution equation 
of the magnetic energy turns out to depend on the cor- 
relation function of the turbulent velocity field. There- 
fore, we present a general model for this correlation func- 
tion for turbulence of arbitrary type with v{£) cx £ in 
the inertial range. In Sec. |III| we calculate the critical 
magnetic Reynolds number for the small-scale dynamo, 
Rm cr i t . Furthermore, we determine the growth rate T 
of the magnetic energy. The latter is calculated in the 
limit Pm — > oo and for finite magnetic Prandtl num- 
bers (Pm > 10). In Sec. IV we discuss our results and 



compare them to the literature. We close our work with 
conclusions in Sec. W\ 



II. ANALYTIC FRAMEWORK 

A. Kazantsev theory 

In 1967 Kazantsev set up a theory for describing the 
time evolution of the magnetic energy that grows due to 
turbulent motions of a conducting fluid |18j . The mecha- 
nism of converting kinetic energy into magnetic energy in 
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this way is known as the turbulent or small-scale dynamo. 
In this section we describe the Kazantsev theory follow- 
ing mainly the formalism proposed by Brandenburg and 
Subramanian [2D] and Subramanian |19j . 



1. The turbulent velocity field 

The theoretical description of turbulence starts with 
decomposition of the velocity field v into a mean field 
(v) and a turbulent component 8v: 

v = (v) + Sv. (1) 

Following the work of Taylor [3S] we model the spatial 
appearance of turbulence via the two-point correlation 
function. The correlation of two turbulent velocity com- 
ponents at the positions ri and r 2 at the times t and s for 
a Gaussian random velocity field with zero mean, which 
is isotropic, homogeneous and (5-correlated in time, is 



(SviiTutyvjfas)) = T tf (r)<f(t - s) 



(2) 



with the two-point correlation function Ty (r), where 
r = |ri — r 2 |. It was shown by Batchelor [35] that the 
correlation function can be divided into a transverse part 
Tn and a longitudinal part XL in the following way: 



Ti 



-W=(^-^)T N (r) + ^T L (r). (3) 



We neglect here the effect of helicity, which would appear 
as an additional term in Ty. 

In the special case of a divergence-free turbulent veloc- 
ity field (div 8v = 0), characteristic of incompressible 
fluids, we find that the transverse correlation function is 
connected to the longitudinal one by 



^(H = ^(r 2 T L (r)) 



(4) 



2r dr 

For the other extreme case, an irrotational turbulent ve- 
locity field (rot Sv = 0), as expected for purely shock- 
dominated flows, we find the relation 



dT N (r) 

Th(r) = r — h X N (r). 

dr 



(5) 



with the two-point correlation function 



M w (r, t) = (<% - r -^pj M N (r, t) + ^Af L (r, *)• (8) 



We will omit the dependencies for a better overview in 
most of the following equations. 

As the magnetic field is always divergence-free, i.e., 
d/dru Mij(r,t) = d/drij Mij(r,t) — 0, we can derive 
a relation between the transverse and the longitudinal 
correlation function similar to Q: 



(9) 



Ml and 



where we have used that (r^j /r 2 )M, 
(ri/rj)Mij = Mn- The time derivative of (SBiSBj) is 

d 



H ((BiBj) - (Bi) (Bj)) 
IdB 



R dB j 



(10) 



In the upper equation we can substitute the induction 
equation 



dB 

~dt 



VxvxB-!)VxVxB, 



(11) 



where 77 = c 2 /(4ira) is the magnetic diffusivity with the 
speed of light c and the electrical conductivity tr, and the 
evolution equation of the magnetic mean field 



9(B) 
dt 



V x [(v) x (B) - ?7 off V x (B) 



with the effective parameter r] c s = rj + 7l(0). 
lengthy derivation (20] this leads to 



(12) 
After a 



9M L A r», 9 ^ 4K diff 
—q^- = 2KdiffM L + 2 



Mi 



Tt. 



Mr 



(13) 



2. Kazantsev equation 

Like the velocity field, the magnetic field can be sepa- 
rated into a mean field (B) and a fluctuation part <5B: 



B = (B) 



SB. 



(6) 



Now let us assume that the fluctuating component SB, 
like the velocity field, is a homogeneous, isotropic Gaus- 
sian random field with zero mean. Then we can write the 
correlation function as 



Mij{r,t) 



(7) 



with 



Kdiff(r) =J7 + T L (0)-T L (r). 



(14) 



The prime denotes differentiation with respect to r. The 
diffusion of the magnetic correlations, Kdiff, contains in 
addition to the magnetic diffusivity r\ the scale-dependent 
turbulent diffusion T L (0) - T L (r). 

With the solution of Eq. (13) we can calculate Mn also 
by using the relation ^ and so find the total correlation 
function of the magnetic field fluctuations My. We note 
that this quantity is proportional to the energy density 
of the magnetic field, B 2 /(8ir). 
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In order to separate the time from the spatial coordinates 
we use the ansatz 



M L (r,i) 



2FI 



Substitution of this ansatz in Eq. (13) gives us 



Kdiff(^) 



dV(r) 



+ U(r)ip(r) = -Tip(r). 



(15) 



(16) 



This is the Kazantsev equation. It formally looks like 
the quantum-mechanical Schrodinger equation with a 
"mass" h 2 /(2Kdis) and the "potential" 

Mr) = <il ^ 2 | 2Kdiff | 2r N , 2 ( y L - T N + « di ff) _ 

(17) 

It describes the kinematic limit, because U is indepen- 
dent of the time. 



B. WKB approximation 

We can use common methods from quantum mechan- 
ics, like the WKB approximation, to solve the Kazantsev 
equation (16). WKB stands for Wentzel, Kramers, and 



Brillouin, who developed this method in 1926 [3"T] . 



1. Solution of the Kazantsev equation in the WKB 
approximation 

In order to use the standard WKB method, we have 
to make some substitutions. Let us first introduce a new 
radial coordinate x by defining r = e x . This leads to 



Kdiff(x) d / 1 dtp(x) 



dx \ e x dx 



{r + U(x))ip(x) =0.(18) 



Next we eliminate the first-derivative terms through the 
substitution 



ip(x) 



W2 



6{x), 



to obtain 



d 2 9{x) 
dx 2 



+ p(x)9(x) = 



(19) 



(20) 



>0.oo 



> 



The boundary conditions for ijj(r) and 9{x) are 
il>(r) 

=> 9{x) r ^ ±0 °) 0. (23) 

We can make some predictions about the shape of the 
function 9(x). For very small x (x — > -co), p(x) goes 
to —1/4 < 0, which leads to exponentially growing and 
decaying solutions of 9. In the other limit (x —> oo), 
p(x) — ¥ —T/rj e 2x , we have growing mode solutions only 
for positive T. The boundary conditions require that 
9 needs to grow exponentially for x — > — oo and decay 
exponentially at x — > oo. In order to arrange this, p{x) 
must go through zero, so U(x) needs to become negative 
for some r. From now on we label the roots of U{x) as 
x\ and X2 > X\. As U(r) becomes negative for some r, 
p(r) becomes positive for certain values of r. This means 
that we have oscillatory solutions for x\ < x < X2- The 
condition for the eigenvalues T in this case is [38] 



ri~i\A ' 2n + 1 
\IP\x ')dx = — 7T 



(24) 



for different excitation levels n £ N. In this work we 
concentrate on the lowest mode n = 0. 



2. Validity of the WKB approximation in general 

In order to find the limits in which the WKB method 
leads to valid solutions of the Kazantsev equation, we 
derive the differential equation that is solved exactly by 



9{x) 



P 



1 / 

l7I exp ±i 



\Jp(x')dx' 



(25) 



The second derivative of 9{x) with respect to x can be 
written as 

3 {pT 



\x) 



1 



P 



Ap 2 16 p 3 



P 9{x) = 0, 



(26) 



where now the prime denotes d/dx. This equation results 
in the Kazantsev equation ( 20 ) if 



1/0*01 « i. 



with 



Ap 2 16 p 3 



l\2 



(27) 



(28) 



with the definition 



We use this result in the Appendix to check the range of 
parameters in which the WKB method produces accurate 
solutions of the Kazantsev equation. 



p(x) 



[r + U{x)]e 2x 



(21) 



The WKB solutions of Eq. ( 20 ) are linear combinations 
of 

9{x) = ^exp Ui J ^pj¥)dx^j . (22) 



C. Modeling the turbulent correlation function 

We analyze the case of general types of turbulence, 
which can be described by the relation between the ve- 
locity v(£) and the size I of a turbulent fluctuation, 



v(£) oc f>. 



(29) 



5 



The power-law index i? varies for the different types, ft 
attains its minimal value of ■& — 1/3 for Kolmogorov 
theory [37], i.e., incompressible turbulence. For Burgers 
turbulence [35], i.e., highly compressible turbulence, i9 
gets its maximal value of 1/2 [55] . 

Motivated by the definition of the scale-dependent tur- 
bulent diffusion coefficient in the last section, 



r?turb(r) =T L (0)-r L (r), 



(30) 



we construct a model for the longitudinal correlation 
function of the turbulent velocity field T^(r). The dif- 



fusion coefficient is calculated from the power law ( 29 ) in 
the following way: 



Tk vab (r)<xv t £oci°t = i° +1 



(31) 



So we assume the correlation function in the inertial 
range to be [191 130] 



r L (r) 



VL 

IT 



(l - (r/Lf^ 



(32) 



The pre factor VL fixes the units, which should be the 
same as for diffusivity. V and L are the velocity and the 
length scale of the largest eddies. On the diffusive scale 
the correlation function should be steadily continued and 
satisfy the condition that its derivative T L (0) vanishes 
at r = 0. This is accomplished, for example, for Tl oc 
r 2 . The exact form of Tl in the diffusive range does not 
affect the results crucially |19j . Furthermore, we expect 
no correlation on scales larger than the largest eddies, so 
Tl should vanish there. 

Taken all together, we can set up a general turbulence 
model for the longitudinal correlation function on the 
different length scales as follows: 



r L (r) 



VL 
3 

VL 
3 



1 



1 - 



Re (i-tf)/(i+tf) (j) 2 ) o < r < 4 

, v (r)" +1 ) U<r<L 
L<r, 

(33) 

where £ c — L Re denotes the cutoff scale of the 

turbulence and L the length of the largest eddies. The 
hydrodynamic Reynolds number Re is defined as (VL)/v 
with the typical velocity of the largest eddies V and the 
viscosity of the gas v. 

The transverse correlation functions Tn for a divergence- 
free (i.e., Kolmogorov turbulence) and for an irrotational 
(i.e., Burgers turbulence) turbulent velocity field can be 
derived from the relations Q and Notice, however, 
that a turbulent velocity field that is divergence free or 
irrotational in the inertial range does not have to be this 
in the diffusive range. We find for the extreme cases in 
the inertial range (£ c < r < L) 



In (r) = 



VL 
~3~ 
VL 
~3~ 



1 - 



1 - 



(z) 
(z) 



4/3 



3/2 



(34) 
(35) 



In order to find a general expression for Tn we make the 
ansatz 



0+1 



(36) 



where t(i9) = a-bd. With Eqs. Q and ([35]) we find that 
a = 21/5 and b = 38/5. Furthermore, we find the small- 
scale transverse correlation (i.e., < r < £ c ) by steady 
continuation. So we end up with the following model for 
the transverse correlation function for the general slope 
of the turbulent velocity spectrum: 



Mr) 



VL 
3 

VL 
3 



1 - iO^Re' 1 -*^ 1 



" ! (zrl ()< r < 4 

£ c < r < L 

i < r, 

(37) 

with = (21 - 38i?)/5. 

The longitudinal and transverse correlation functions de- 
pend on the dimensionless parameter y = r/L as shown 
in Fig. [TJ for Kolmogorov and Burgers turbulence. We 
choose here a fixed hydrodynamical Reynolds number of 
10 5 . In the inset of Fig. [l] we show a zoom into the dis- 
sipative range (0 < r < £ c ). Furthermore, in Fig. [2] we 
plot the potential of the Kazantsev equation, resulting 
from our model of Tl for Kolmogorov and Burgers tur- 
bulence. We choose two different values for the Reynolds 
number Re and different magnetic Prandtl numbers Pm. 
The magnetic Prandtl number is defined as the ratio of 
the magnetic Reynolds number Rm = (VL)/rj and the 
hydrodynamical Reynolds number Re. The potential at 
fixed Re and Pm is deeper in the small-scale range in the 
Kolmogorov case than in the Burgers case. For higher 
Re the potential gets deeper and the cutoff scale £ c de- 
creases. For higher Pm the potential in the small-scale 
range gets broader. 



III. SOLUTION OF THE KAZANTSEV 
EQUATION 

In this section we use our model of the turbulent veloc- 



ity correlation function ( 33 ) and ( 37 ) as the input for the 
Kazantsev theory. We solve the Kazantsev equation in 
order to obtain the characteristic properties of the small- 
scale dynamo. We use the WKB method, which gives 
a good approximate solution for large magnetic Prandtl 
number. In fact, in the limit of infinite magnetic Prandtl 
number the WKB approximation is an exact solution of 
the Kazantsev equation. For more details about the va- 
lidity of this approximation, see the Appendix. 



A. Critical magnetic Reynolds number for 
small-scale dynamo action 

Intuitively, one expects that the high magnetic diffu- 
sivity for very low magnetic Reynolds numbers prevents 
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FIG. 1. (Color online) Dependence of the longitudinal and 
transverse correlation functions TL and Tn on the dimension- 
less parameter y = r/L for Kolmogorov (1? = 1/3) and Burg- 
ers ($ = 1/2) turbulence. We choose a fixed Reynolds number 
of 10 s . The vertical lines indicate the cutoff scale of the tur- 
bulence £ c and the largest scale of the eddies L. Notice that 
the cutoff scale for Kolmogorov turbulence (£^ = Re _3//4 L) is 
different from the one for Burgers turbulence (£™ = Re~ 2/3 L). 
The inset shows a zoom of the dissipative range. 



amplification of the magnetic field. Even higher diffu- 
sivity eventually results in a net decrease of the field 
strength. In this section we calculate the critical mag- 
netic Reynolds number Rm cr i t for small-scale dynamo 
action. To accomplish this we set the growth in our equa- 
tions at zero. 

It should be noted that we use the inertial range (£ c < 
r < L) for determining Rm cr j t as the potential in this 
range is always negative and for that we have a positive 
growth rate (see Fig. [2| . In this range with our turbu- 
lence spectrum and T — we get for the p- function (21) 



p(y) 



-9/4 - a (tf)Rm cl ,t/ +1 + b(tf)Rmi it yW+V 



[1 + \Bm cll tV d+1 ) 



(38) 



with a(tf) = 5/6 - (79/30) + (157/30) 1? 2 and 6(i?) 
(14/15) d- (103/60) z9 2 . 



Now we can evaluate the eigenvalue condition (24) for 
this p(y) in the "ground state" n = 0: 



Vp(y)— = 
wo, y z 



(39) 



in which the additional y comes from the substitution 
y = r/L = e x /L. The limits of the integral are the roots 
of p{y). There is only one real and positive root of p(y), 
which we label y±. For the upper limit we have to realize 
that the potential (17) changes for y > 1 to 2i]/(yL) 2 , 
which is clearly always positive. Furthermore, also the 
diffusion coefficient K di g = 77 + T L (Q) > for y > 1. With 
U and Kdiff being positive p{y) is negative in this range, 
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FIG. 2. (Color online) Dependence of the potential on the 
dimensionless parameter y = r/L for Kolmogorov (i? = 1/3) 
and Burgers (i9 = 1/2) turbulence. We choose two differ- 
ent Reynolds numbers Re = 10 5 and Re = 10 7 , and differ- 
ent Prandtl numbers Pm = 10 2 , Pm = 10 4 and Pm = 10 8 . 
The cutoff scale l c depends on the turbulence model and the 
Reynolds number. For Kolmogorov turbulence l c = Re _3//4 L; 
for Burgers turbulence £ c = Re~ 2 / 3 L. A Reynolds number 
10 x is indicated in the cutoff scale as 



which means that p{y) needs to go through zero during 
this transition. So we have our second root at roughly 
r « L and y% = 1. 



We can solve Eq. (39) numerically for the critical mag- 
netic Reynolds number Rm cr i t if we put in a fixed value 
of i?. Recall that d was defined in the inertial range 
of the turbulence via the relation v(£) oc £ . Results 
for common models in the literature can be found in 
Table |TJ In Fig. [3] we show how the critical magnetic 
Reynolds number depends on Here one can see that 
the critical magnetic Reynolds number increases rapidly 
as gets closer to its maximum value of 1/2. An em- 
pirical fit Rm cr it,fit($) through these data in the range 
0.33 < < 0.5 is 



Rmcrit,fit(tf) = 88 (tan(2.68tf + 0.2) - 1) . 



(40) 



Furthermore, we collect the results for common turbu- 
lence models in the literature in Table [T| 
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We find that the small-scale dynamo is more easily ex- 
cited in the case of a purely rotational turbulent veloc- 
ity field, i.e., for Kolmogorov turbulence, where we find 
Rm cr it 110. The critical magnetic Reynolds number 
for a turbulent field with a vanishing rotational compo- 
nent, i.e., Burgers turbulence, is roughly 2700. 
The results are discussed in Sec. HVl 
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FIG. 3. (Color online) Dependence of the critical magnetic 
Reynolds number Rm cr i t from the slope of the turbulent ve- 
locity spectrum i?. The dashed line is an empirical fit through 
our results. 



B. Growth rate of the small-scale magnetic field 

1. Growth Rate in the Limit Pm — ?> oo 

In this section we derive a general analytical solution 
for the growth rate T for an arbitrary slope of the tur- 
bulent velocity spectrum, in the limit of infinite Prandtl 
number. 

As the potential has its minimum in the small-scale 
range, i.e., the dissipative range of the turbulence (see 
Fig. [2]), the growth rate, which is the eigenvalue of the 
Kazantsev equation, takes its maximal value there. So 
we expect the fastest growing mode to be in the small- 
scale range. 

In order to have scale-independent equations, we intro- 
duce the substitution 



WRe 
3L V 



1/2 



'Re 3/2 Pm~ 



1/2 



(41) 



where the magnetic Prandtl number is Pm=Rm/Re. 
The p function in z space [see (21 ) and (41 )] for the gen- 



eral turbulence spectrum ( 33 ) in the dissipative range is 
p(z) = 



A z* - B z 2 - 45Re( 3+7 "^ 2+2 "> 



20Re 



x ' 2 (r^ 1 " 



■3i5)/(2+2i?) 



Re 



1/(1+*) ^2 



(42) 



with the definitions 



20, 



A = Re (5+")/(2+2tf) (163 _ 30 4tf) _ ^Re 5 / 2 f , (43) 
B = (304tf - 98) Re 2 + ^Re< 2+M >/< 1+ *>r , (44) 



and the normalized growth rate 

V 



(45) 



In the limit of large Prandtl number z is large, too, and 
we can neglect the constant terms. We obtain 



p{z) 



Re -(5+*)/(2+2tf) Aq . 



B 



20 



(46) 



The one real and positive root of this function is Z\ — 
s/Bq/Aq. At the cutoff scale of the turbulence the p 
function changes its sign. We take this as our second 
root and so have z 2 = R e (3, '- 3)/(4 ' ?+4) . So we 

get for the general eigenvalue condition 



Re -(5+tf)/(4+«) fZ2 I Aqz2 _ Bq ^ 

27! L V^^ dz =2< 

resulting in the analytical solution of the integral 

R e -(5+i?)/(4+4i?) 

2V5z 



A ln ( 2 ( \fA z+ 

Z2 7T 



(47) 



y/A z 2 -Bo)) -y/A z 2 -B 



(48) 



For ^2^1 this becomes 
Re -(5+tf)/(4+4#) 

2~vi A 



A 



1 - In U^fAo 



)Z2 



+ -ln(4B ) 



7T 

2' 



(49) 



A zero-order iterative solution for T gives us 
163- 304i?, 



r = 



60 



-Re 



^ Re -(5+*)/(4+4^) ^ _ ln (4^^) + 1/2 In (4B )] 
which becomes for large Prandtl number 
- 163 - 304$ 



60 



-Re (1 ~' 5)/(1+ ' 5) . 



(51) 



As a result we get for the absolute growth rate V for a 
general slope of the turbulent velocity spectrum 

r- (163 ~y ^ Rc (i-*)/a+*) (52) 
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FIG. 4. (Color online) The normalized growth rate of the 
small-scale dynamo in the limit of infinite magnetic Prandtl 
number, depending on the Reynolds number Re. For the 
slopes of the turbulent velocity spectrum i? we choose com- 
mon values from the literature: K41 [27|, SL94 [29], BNP02 
[H], L81 [3Tj, FRKSM10 [32] (sol: solenoidal forcing; comp: 
compressive forcing), OM02 [33] and B48 [28) . 



FIG. 5. (Color online) Dependence of the normalized growth 
rate of the small-scale dynamo on the magnetic Prandtl num- 
ber Pm for Kolmogorov turbulence. We choose different val- 
ues of the Reynolds number. Notice that in the limit Pm — > oo 
the normalized growth rates are f = 10.28 for Re = 10 2 , 
f = 102.78 for Re = 10 4 , and f = 10277.78 for Re = 10*. 
These limits are indicated in the plot as horizontal lines. 



in the limit Pm — > oo. In Fig. [3] we show the dependency 
of the normalized growth rate T on the Reynolds num- 
ber for different types of turbulence. One extreme case is 
incompressible turbulence, i.e., Kolmogorov turbulence, 
with r oc Re ' . In the other extreme case, highly com- 
pressible turbulence, i.e., Burgers turbulence, the growth 
rate increases only as Re 1 / 3 . Altogether we find that the 
growth rate increases faster with the Reynolds number 
when the compressibility is lower. 



2. Growth rate as a function of the Prandtl number 

In this section we discard the assumption of infinite 
Prandtl number. In this case we have to solve the full 



equation resulting from the WKB method ( 24 ) 



z a - 2 : 



(53) 



with p(z) from (42 1. There is no analytical solution of 
this integral equation. 

The numerical results of the normalized growth rate are 
shown in Fig. [5] for Kolmogorov turbulence and in Fig. [6] 
for Burgers turbulence. We plot the normalized growth 
rate depending on the Prandtl number for different values 
of the Reynolds number. 



IV. 



DISCUSSION 



The main results of this work are the critical magnetic 
Reynolds number Rm cr j t and the growth rate T of the 
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FIG. 6. (Color online) Dependence of the normalized growth 
rate of the small-scale dynamo on the magnetic Prandtl num- 
ber Pm for Burgers turbulence. We choose different values of 
the Reynolds number. Notice that in the limit Pm — > oo the 
normalized growth rates are T = 0.85 for Re = 10 2 , T = 3.95 
for Re = 10 4 , and f = 85.1 for Re = 10 s . These limits are 
indicated in the plot as horizontal lines. 



small-scale dynamo given by Fig. H and (52}. In Table [I) 
we collect our results for typical turbulence models from 
the literature, defined via the relation v(£) oc £ . The two 
cxtrema of turbulence are Kolmogorov turbulence [27] . 
i.e., incompressible turbulence, with i? = 1/3 and Burg- 
ers turbulence with $ = 1/2, which describes highly com- 
pressible turbulence with vanishing rotational component 
[28] . Between those extreme values we choose values of 
d from observations of molecular clouds, like d w 0.38 
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from Larson [3T] and i? sa 0.47 from Ossenkopf and Mac 
Low |33j . Furthermore, we give $ ~ 0.35 as an example 
for a theoretical model of intermittency [53] . Numerical 
experiments give $ « 0.37 for driven supersonic magne- 
tohydrodynamical (MHD) turbulence [30] and z9 « 0.43 
and $ « 0.47 for solenoidal and compressive forcing of 
the turbulence 32J. Notice, however that the mean val- 
ues of $ from observations and simulations have a typical 
uncertainty of 10%. 



A. Critical magnetic Reynolds number 



We have evaluated Eq. ( 39 ) for different slopes of the 



turbulent velocity spectrum The results for common 
types of turbulence from the literature are listed in Table 
[I] and pictured in Fig. [3] 

From these results we see that for all types of turbulence 
a high magnetic Reynolds number needs to be exceeded 
for small-scale dynamo action. In astrophysical objects 
we often find very high magnetic Reynolds numbers (see 
the compilation in |41j). The core of Jupiter, for example, 
has Rm w 10 6 , the solar convection zone has Rm ss 10 s , 
and the solar corona already Rm sa 10 12 . In the interstel- 
lar medium we find Rm 10 17 and in a typical galaxy 
about 10 19 . Consequently, the critical magnetic Reynolds 
number is exceeded by far in nature, and we expect that 
a small-scale dynamo is excited in typical astrophysical 
objects. 

There are different ways to obtain approximate solutions 
of the Kazantsev equation. In addition to the WKB 
method there is also an asymptotic solution, which uses 
the separation of scales [24] • The potential U(r) and the 
mass m(r) are estimated on different scales, such that a 
solution of the Kazantsev equation can be found. Ro- 
gachevskii and Kleeorin [22] used this method to deter- 
mine the critical magnetic Reynolds number for differ- 
ent compressibilities of turbulence. They found that, in 
the limit of small magnetic Prandtl number Rm needs 
to be larger than roughly 400 for excitation of a small- 
scale dynamo in the case of Kolmogorov turbulence. For 
a larger compressibility, i.e., toward Burgers turbulence, 
they found that Rm cr i t increases sharply. We see the 
same trend of increasing Rm cr i t for higher compressibil- 
ity. However, our result for Kolomogorov turbulence dif- 
fers by a factor of roughly 4, which might be caused 
by the fact that we analyze the limit of large magnetic 
Prandtl number. 

Recent high-resolution numerical studies confirm the ex- 
istence of a critical magnetic Reynolds number for small- 
scale dynamo action. Haugen et al. [42] found Rm cr it ~ 
35 for subsonic turbulence and Rm cr i t ~ 70 for super- 
sonic turbulence at a magnetic Prandtl number of about 
unity. In numerical simulations, the magnetic Reynolds 
number can be estimated by Rm sa (A/£ c ) +1 , where A is 
the typical size of turbulent structures and £ c is the cut- 
off scale of the turbulence. The latter can be estimated 
by l c w 0.5 Ax with Ax the minimal resolved size in a 



simulation 43J. In resolution studies, Sur et al. [T2] and 
Federrath et al. [13] found that the typical length of a 
turbulent fluctuation needs to be resolved with at least 
30 grid cells in magnetohydrodynamic simulations of a 
self- gravitating gas. Only then is the magnetic field am- 
plified exponentially, which is explained by the action of 
the small-scale dynamo. 

For a physical interpretation of this result it is useful to 
take the stretch-twist- fold- dynamo as a toy-model of the 
turbulent dynamo. In this picture we think of a magnetic 
flux rope that gets stretched due to turbulent motions. 
This motion decreases the cross section of the flux rope. 
If we assume that the magnetic flux is frozen in the gas, 
then the magnetic field increases during this process, be- 
cause the magnetic flux is a conserved quantity [2H1 SI] • 
This process works best in a purely rotational turbulent 
velocity field. Therefore, we expect the dynamo to be 
more easily excited in Kolmogorov turbulence. In order 
to see this process in simulations, one needs to resolve the 
stretching, twisting, and folding of the field lines, which 
explains the required high resolution. 
The determination of the critical magnetic Reynolds 
number is, moreover, the first step to understanding the 
saturation of the small-scale dynamo. For if the magnetic 
field in a system increases, back reactions from the gas 
become more important. Then processes like ambipo- 
lar diffusion can change the properties of the gas and the 
magnetic Reynolds number can decrease. If the magnetic 
Reynolds number becomes smaller than the critical mag- 
netic Reynolds number, the magnetic field stops growing 
and the small-scale dynamo is saturated. 



B. Growth rate 

In Fig. [3] as well as in Table [I] we present our results 
for the growth rate of the small-scale dynamo in the limit 
of infinite magnetic Prandtl number. Our results show 
that the growth rate is proportional to the velocity V of 
the largest eddy divided by its length L. The ratio V/L 
is the reciprovcal of the turnover time of an eddy. Thus, 
the growth rate increases with decreasing turnover time, 
and the smallest modes grow at the highest rate. This is 
expected, because smaller turnover times lead to a faster 
tangling of the magnetic field lines. 

Furthermore, the growth rate increases with increasing 
hydrodynamical Reynolds number for all types of turbu- 
lence, characterized by v(£) oc I . In order to achieve 
the same growth rate for Kolmogorov and Burgers tur- 
bulence we have to provide a larger Reynolds number 
in the latter case. Assuming a fixed Reynolds numbers 
in both extreme cases, Re and Re , the growth rates 
of the two different turbulence types are the same for 
Re F 



0.18(Re B ) 3 / 2 . This fact can again be explained 
with the stretch-twist-fold model (see Sec. IV A). We 
need solenoidal modes, i.e., divergence- free modes, of 
the turbulence for this process |34j . which explains why 
incompressible turbulence amplifies the magnetic field 
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Model and reference 


c) 


Rm crit 


r (Pm -> oo 


Kolmogorov [27] 


1/3 


« 107 


I Re 1/2 

36 


Intermittencv of Kolmogorov turbulence (She and Leveque 29 ) 


0.35 


« 118 


0.94 Re 48 


Driven supersonic MHD-turbulence (Boldyrev et al. 30 ) 


0.37 


« 137 


0.84 Re 46 


Observation in molecular clouds (Larson [31]) 


0.38 


« 149 


0.79 Re 


Solenoidal forcing of the turbulence (Federrath et al. |32|) 


0.43 


« 227 


0.54 Re 40 


Compressive forcing of the turbulence (Federrath et al. [32]), 


0.47 


« 697 


0.34 Re 36 


Observations in molecular clouds (Ossenkopf and Mac Low [33]) 









Burgers [25] 1/2 ^2718 |i Re 1/3 



TABLE I. The critical magnetic Reynolds number Rm cr i t and the normalized growth rate of the small-scale dynamo F in the 
limit of infinite magnetic Prandtl number. We show our results for different types of turbulence, which are characterized by 
the exponent -d of the slope of the turbulent velocity spectrum, v(i) oc £ . The extreme values of •& are 1/3 for Kolmogorov 
turbulence and 1/2 for Burgers turbulence. 



more effectively. 

With the asymptotic solution of the Kazantsev equation, 
Rogachevskii and Kleeorin [22] found in the limit of small 
magnetic Prandtl number T oc ln(Rm/Rm cr ;t). The con- 
stant of proportionality depends on the amount of com- 
pressibility. In a later work Kleeorin and Rogachevskii 
[2"5] found that this logarithmic scaling of the growth 
rate is valid only in the vicinity of the threshold of small- 
scale dynamo excitation. For magnetic Reynolds num- 
bers much larger than Rm cr i t , they found in the limit 

1 /9 

of small magnetic Prandtl number T oc Rm ' for Kol- 
mogorov turbulence. As for a constant magnetic Prandtl 
number Rm oc Re, this agrees with our result. 
There are recent high-resolution numerical simulations 
that model the turbulent dynamo. The two limiting ways 
of driving the turbulence are solenoidal, i.e., divergence- 
free, forcing and compressive, i.e., rotation-free, forc- 
ing. These simulations show, in agreement with our 
study, that solenoidally driven turbulence leads to larger 
growth rates of the small-scale dynamo. Waagan et al. 
[4"4"j found, using a Reynolds number of about 1500 [4"5] . 
and a magnetic Prandtl number of about 1, for totally 
solenoidal forcing of the turbulence f so1 = 0.60 and for 
totally compressive forcing r comp — 0.28. These values of 
the growth rate are about a factor of 17 lower than those 
from our model (with Re = 1500), f so1 = f^ 0A3 10.07 
and r comp = r^ - 47 w 4.73. This can be explained by 
the fact that the simulations had a very low magnetic 
Prandtl number of about 1. However, our result for the 
growth rate in Table [I] was derived with the assumption 
of infinite Prandtl number. We have also explored the 
range of lower Prandtl numbers. The result is presented 
in Fig. [5] for Kolmogorov turbulence and in Fig. [6] for 
Burgers turbulence. But with our model, we can make 
no predictions for Prandtl numbers around unity, because 
in this range the WKB approximation is no longer ap- 



plicable (see the Appendix). However, the trend is that 
the growth rate decreases for lower Prandtl numbers and 
this can explain the lower growth rates from the simu- 
lations. Yet the ratio of the growth rate of turbulence 
driven by solenoidal and compressive forcing is in both 
cases about 2_(our model: f sol /f comp « 2.1; Waagan et 
al. [H]: r sol /r comp ~ 2.1), which shows that incompress- 
ible turbulence is more efficient in amplifying a magnetic 
field via the small-scale dynamo. 

Furthermore, Federrath et al. [34] have presented a study 
of the Mach number dependency of the growth rate of 
the small-scale dynamo, where they compared solenoidal 
with compressive forcing of the turbulence. They found 
that, for low Mach numbers, the ratio of the growth rate 
of turbulence driven by solenoidal and compressive forc- 
ing is about 30 (for Mach number M = 0.1: f so1 rj 1.2, 
pcomp ^ 0.04). However, for higher Mach numbers their 
calculations also result in a ratio of f sol /f comp w 2 (for 
Mach number M = 10: f so1 « 0.7, f comp « 0.3), which 
is in agreement with our results. The lower growth rates 
in the simulation again may come from low Prandtl num- 
bers in the simulations, of the order of unity. However, 
a great uncertainty is the Reynolds number in numerical 
simulations, which is only a crude estimate. 



V. SUMMARY AND CONCLUSIONS 

We have presented an analytical treatment of the 
small-scale dynamo, using the Kazantsev theory. For 
that we have modeled the correlation function of the 
turbulent velocity field, depending on the slope of the 
turbulent velocity spectrum, fl, in v(£) oc £ . With this 
model, we solved the Kazantsev equation in the WKB 
approximation and tested the validity of this approxi- 
mation. We determined the critical magnetic Reynolds 
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number for the small-scale dynamo and its growth rate in 
the case of infinite and finite magnetic Prandtl numbers. 
The main results of our work are as follows: 

• The critical magnetic Reynolds number Rm cr i t for 
the small-scale dynamo increases as the exponent i? 
increases (see Fig. [3]). For Kolmogorov turbulence 



Rm- 
2700. 



110 and for Burgers turbulence Rm c 



The growth rate of the magnetic field energy in the 
limit of infinite magnetic Prandtl number is 



r 



(163 - 304^) V Re(1 ^ )/(1+1?) 



60 L 

(see also Figure [4]). 



(54) 



• For decreasing magnetic Prandtl number the 
growth rate decreases. The details of this drop de- 
pend on the type of turbulence (see Figs. [5] and 

©■ 

• A validity test shows that the WKB-approximation 
gives exact solutions in the limit of infinite mag- 
netic Prandtl number. The approximation breaks 
down at a Prandtl number of around unity (see 
Figs. |7j and ||. 

With these results we are able to make predictions about 
the first magnetic fields in the Universe. As the hydro- 
dynamical Reynolds numbers are typically very high in 
astrophysical objects, the growth rates of the small-scale 
dynamo are very high, too. We thus expect that the 
time until saturation is shorter than the collapse time of 
a halo 13 . For that reason we might already have high 
magnetic field strengths even before the formation of the 
first stars, the first galaxies, and the first galaxy clusters. 
Turbulence and magnetic fields are key ingredients of 
current star formation theory [JSJ H?]. Magnetic fields 
drive jets and outflows from young stars. Stellar winds 
and supernova explosions, which end the lives of mas- 
sive stars, enrich the interstellar medium with heavy el- 
ements forged in the stellar interior. These processes are 
crucial for the chemical composition of the Universe, de- 
termining cooling and heating processes in the gas. This, 
in turn, is very important for the formation of the next 
generation of stars. The momentum from jets and out- 
flows around accreting protostars may disperse some of 
the envelope material that otherwise would fall onto the 
central star. Thus, they are important ingredients for 
our understanding of the physical origin of the observed 
distribution of stellar masses |48] . 



ACKNOWLEDGMENTS 

J.S. acknowledges financial support by the Deutsche 
Forschungsgemeinschaft (DFG) in the Schwerpunktpro- 
gramm SPP 1573 "Physics of the Interstellar Medium" 



under Grant No. KL 1358/14-1. D.R.G.S. acknowl- 
edges funding through the SPP 1573 (Project No. 
SCHL 1964/1-1) and the SFB 963 "Astrophysical 
Flow Instabilities and Turbulence". C.F. acknowl- 
edges funding from the Australian Research Council 
(Grant No. DP110102191) and from the European Re- 
search Council (FP7 Grant Agreement No. 247060). 
C.F., R.B., and R.S.K. acknowledge subsidies from the 
Baden- Wiirttemberg-Stiftung under Research Contract 
No. P-LS-SPII/18 and from the German Bundesminis- 
terium fur Bildung und Forschung via the ASTRONET 
project STAR FORMAT (Grant No. 05A09VHA). R.S.K 
also thanks the DFG for financial support via Grants 
No. KL1358/10 and No. KL1358/11, as well as via the 
SFB 881 "The Milky Way System". R.B. acknowl- 
edges funding by the Emmy-Noether Grant No. (DFG) 
BA 3706 and the DFG via the Grands No. BA 3706/1-1 
and No. BA 3706/3-1. 



APPENDIX: VALIDITY OF THE WKB 
APPROXIMATION 

The WKB method is only an approximate solution of 
the Kazantsev equation. We have derived condition (27 1, 
|/| <C 1, for which the WKB-method is valid in order to 
find solutions. In z space, / reads 



/(*) 



z 2 p"{z) + 2zp'{z) 3 [zp'{z)f 



Ap{z) 



16 p(z) 3 



(55) 



However, we have seen that the magnetic field is ampli- 
fied most strongly on the scale £ c (z) = ^/Pm/3, as here 
the potential U has its minimum. So we analyze f(z,T) 
on this scale and get a dependency on the Prandtl num- 
ber Pm. Hence we label /(4,T) = /(Pm,r). 
One can show that /(Pm, V) vanishes in the limit of large 
Prandtl number for all T and all turbulence types, 



lim /(Pm,r)=0. 

Pm— 7-oc 



(56) 



This means that the WKB method is very good in the 
limit of large magnetic Prandtl number. 



1. Validity of the WKB approximation for 
Kolmogorov turbulence 

In order to check also lower Prandtl numbers we plot 
/(Pm, r) for different normalized growth rates T (45| and 
Kolmogorov turbulence in Fig. [7] However, one can show 
that /(Pm, T) does not depend on the Reynolds number 
for Kolmogorov turbulence. So we choose values for f be- 
tween and the maximal value f max for the plot in Fig. [7J 
where r max is the value for an infinite Prandtl number 
and depends on the Reynolds number. One can see that 
the critical Prandtl number for the WKB approximation 
gets larger with increasing normalized growth rate. 
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FIG. 7. (Color online) The function /(Pm,f) for different 
values of the normalized growth rate for Kolmogorov turbu- 
lence. r max is the normalized growth rate in the limit of in- 
finite magnetic Prandtl numbers, F max = (37/36)Re 1//2 (see 
Sec. Ill for the derivation). The WKB approximation is valid 
within the non-hatched area, i.e., for |/(Pm, T)\ < 0.1. 



To make a more quantitative estimate of the criti- 
cal Prandtl number, we have hatched the area above 
/(Pm,f) = 0.1 and below /(Pm,f) = -0.1. When / 
is not in this area its absolute value is smaller than 10% 
of 1. We take this as a threshold for our approximation. 
We find that our method is applicable in the case of F = 
for 



Pm > 13. 



(57) 



For higher normalized growth rates the critical Prandtl 
number increases. 



Validity of the WKB approximation for Burgers 
turbulence 



We can analyze the validity of the WKB solutions for 
Burgers turbulence in the same way as for Kolmogorov 
turbulence using criterion (27). 

However, we find that the function / given in (|55|) now 



depends not only on the normalized growth rate T and 
the Prandtl number Pm, but also on the Reynolds num- 
ber Re. The result is shown in Fig. |8j where we plot / 
against the Prandtl number for different Reynolds num- 
bers and different normalized growth rates. 
We again determine the critical Prandtl number for the 
WKB method for a vanishing normalized growth rate. 
For our different values of the Reynolds number we 
get the following critical Prandtl numbers at vanishing 



Pm(Re = 10 2 ) > 500, 
Pm(Rc = 10 4 ) > 1100, 
Pm(Rc = 10 8 ) > 5100. 



(58) 
(59) 
(60) 
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FIG. 8. (Color online) The function /(Pm,r) for fixed 
Reynolds numbers and different values of the normalized 
growth rate for Burgers turbulence. Notice that in the limit 
Pm — > oo the normalized growth rates are P = 0.85 for 
Re = 10 2 , F = 3.95 for Re = 10 4 , and P = 85.1 for Re = 10 8 . 
The WKB approximation is valid within the non-hatched 
area, i.e., for /(Pm, V) < 0.1. 
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